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ABSTRACT 


Models  of  human  search  and  target  acquisition  are  typically  based  on  the 
assumption  of  an  independent  glimpse  probability  leading  to  an  exponential  model  of  target 
acquisition  probability.  This  can  be  extended  to  the  case  of  obscurants  which  intermittently 
interrupt  the  line  of  sight  to  any  point  of  the  image  field  as  might  be  the  case  with  drifting 
smokes.  The  model  was  previously  addressed  by  Monte  Carlo  methods  by  Kowalczyk 
and  Rotman  and  here  is  given  a  complete  analytic  solution  by  comparison  with  physical 
models  of  simple  ferromagnets.  The  result  is  that  the  single  exponential  characteristic  of 
the  unobscured  search  model  is  replaced  with  two  exponentials  whose  decay  times  and 
amplitudes  are  determined  by  the  average  duration  times  for  clear  and  obscured  lines  of 
sight. 

The  modified  search  model  still  assumes  that  the  target  acquisition  task  can  be 
decomposed  into  short  glimpses.  If  a  fixed  moderate  length  duration  of  unobscured  time  is 
needed  for  the  completion  of  a  task  (such  as  preparing,  firing,  and  guiding  a  weapon),  a 
different  calculation  must  be  made.  Using  the  obscurant  model  developed  for  the  search 
task,  the  timelines  for  such  situations  are  given  analytically. 
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EXECUTIVE  SUMMARY 
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The  assessment  of  the  effects  of  smoke  on  the  battlefield  presents  a  uniquely 
difficult  problem  for  war  gamers.  In  particular,  simulations  such  as  JANUS  or 
CASTFOREM  play  smoke  at  a  level  of  fidelity  such  that  the  players  affected  by  a  given 
smoke  round  would  be  identified  specifically,  but  the  details  of  the  exact  nature  of  the  effect 
would  be  treated  statisdcally.  In  other  words,  while  these  simulations  are  high  fidelity, 
they  do  not  simulate  a  "virtual  reality"  in  which  the  details  of  the  smoke  plume,  such  as 
fluctuations  and  holes,  are  represented.  The  incorporation  of  such  a  high  level  of  detail  is 
inconsistent  with  the  design  philosophy  of  these  simulation  tools.  The  level  of  hardware 
and  software  reconEguration  that  would  be  required  to  support  such  detail  would  be 
disrupdve  to  say  the  least 

Yet  the  nature  of  the  target  acquisidon  task  must  be  affected  by  these  small  scale 
variadons.  The  existence  of  fluctuadons  and  holes  in  smoke  plumes  has  long  been  known, 
but  has  only  recently  been  explained  quantitadvely.  How  can  these  effects  be  both 
accurately  and  economically  incorporated  into  the  war  games?  This  is  the  quesdon  which 
we  answer  in  this  paper. 

The  tiamewoik  of  our  analysis  is  a  model  of  target  acquisidon  in  smoke  which  was 
developed  previously  by  Rotman.  In  his  treatment,  Rotman  demonstrates  the  appropriate 
averaging  procedure  for  target  acquisition  in  smoke.  The  analysis  herein  accomplishes 
three  things.  First,  we  justify  the  model  by  establishing  the  connecdon  between  it  and 
some  new  results  in  the  physics  of  smoke  plumes.  Second,  we  solve  it  by  deriving  simple 
analyde  soludon  to  the  model.  Third,  we  extend  it  by  showing  how  the  model  applies  to 
the  computation  of  the  probability  of  successfully  completing  a  command  guidance  task. 

The  justification  of  the  original  model  addresses  two  issues.  In  both  cases  we 
show  how  to  determine  the  inputs  required  by  the  Rotman  model  from  the  readily  measured 
properties  of  smoke  plumes.  The  first  issue  is  the  structure  of  the  fluctuations  in  smoke 
plumes.  In  discussing  this  issue  we  rely  on  vory  recent  work  done  by  the  Army  Battlefield 
Environments  Directorate,  which  applies  the  theory  of  atmospheric  turbulence  to  smoke 
plumes.  The  second  issue  is  the  main  simplifying  assumption  of  the  model,  namely  that 
for  purposes  of  target  acquisition,  a  given  line  of  sight  is  either  clear  or  obscured. 
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This  observation  draws  on  our  everyday  experience;  smoke  plumes  seem  to  have  fairly 
well  defined  edges,  and  the  fluctuations  in  concentration  are  evident  to  us  as  actual  holes  in 
the  plume.  We  show  that  this  intuitive  assessment  is  borne  out  by  computations  based 
upon  the  Night  Vision  Static  Model,  which  shows  that  under  many  conditions  (which  we 
determine),  the  predicted  target  acquisition  performance  changes  rapidly  over  a  narrow 
range  of  smoke  concentration. 

The  second  accomplishment,  namely  the  solution  of  the  Rotman  model,  is  based  on 
the  mathematical  similarity  of  this  model  to  a  certain  physical  model  of  ferromagnetism. 
Even  given  the  physical  paradigm,  though,  much  original  work  still  needed  to  be  done 
since  the  relevant  questions  for  the  present  case  do  not  map  onto  the  ones  asked  in  the 
physical  problem.  So  although  the  arithmetic  does  sometimes  gets  thick,  this  complexity  is 
not  gratuitous  because  the  work  does  not  exist  elsewhere.  Fortunately,  the  intermediate 
complexity  can  be  put  aside  at  the  end  since  the  final  result  is  quite  simple. 

Finally,  the  extension  of  the  model  addresses  the  operational  question  that  follows 
naturally  from  the  acquisition  phase:  Given  that  we  have  computed  the  probability  of  target 
acquisition,  can  we  use  the  same  framework  to  compute  the  probability  of  maintaining  an 
uninterrupted  line  of  sight  to  the  target  for  the  time  that  is  needed  to  successfully  engage  it? 
The  answer  is  found  to  be  yes,  but  some  approximations  need  to  be  introduced  if  the 
solution  is  to  be  in  closed  form. 

Further  study  is  needed  before  insertion  into  the  war  game  simulations  should  be 
considered  seriously.  Apart  from  obvious  concerns  over  validation,  the  model  should  be 
exercised  in  a  sensitivity  analysis,  to  determine  whether,  and  under  what  circumstances,  the 
new  predictions  affect  the  determination  of  relevant  measures  of  effectiveness  in  the  war 
games.  This  issue  continues  to  be  addressed  in  the  Army's  Acquisition  and  Simulation 
(AfTQSIM)  Program.  Apart  from  its  possible  insertion  into  the  war  games,  the  model 
developed  by  Rotman  and  improved  here  is  available  as  a  stand-alone  tool  for  back-of-the- 
envelope  assessments. 
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I.  INTRODUCTION 


Rotman^  has  merged  a  simple  model  of  target  acquisition  with  a  simple  model  of 
obscuration.  He  assumes  a  simple  on-or-off  model  of  smoke  obscuration,  and  that  the 
probability  of  target  not  acquired  is  given  by  a  decaying  exponential  in  the  unobscured  time 
on  target.  The  interesting  and  nontrivial  result  is  that  in  order  to  predict  the  mean  time 
required  to  acquire  a  target  in  the  dynamic  environment,  it  is  insufficient  to  simply  specify 
the  level  of  obscuration;  it  is  also  necessary  to  know  something  about  the  temporal 
correlations  of  the  obscured  condition.  In  other  words,  the  e^ect  of  smoke  clouds  that  last, 
say,  2  seconds  with  2-second  intervals  between  clouds  is  much  different  than  those  that  last 
20  seconds  with  20-second  inter  als,  even  though  the  average  transmission  values  are  the 
same  for  the  two  cases. 

Rotman's  model  results  were  obtained  via  Monte  Carlo  simulation  of  the  obscurant 
condition.  We  have  found  that  it  is  possible  to  obtain  an  analytic  solution  of  his  model. 
This  has  the  obvious  benerit  that  it  now  becomes  possible  to  obtain  certain  results  without 
recourse  to  simulations.  An  additional  benefit  that  accrues  may  in  fact  be  more  important: 
The  methodology  which  we  introduce  for  handling  the  smoke  model  is  sufficiently  general 
that  we  can  pose — and  answer — questions  outside  the  realm  of  the  simple  search  model; 
for  example,  questions  related  to  the  effect  of  smoke  on  the  subsequent  engagement  of  the 
target. 

In  this  paper  we  develop  an  analytical  formalism  for  the  solution  of  Rotman's 
model,  and  apply  it  to  two  problems  of  military  interest:  modeling  the  effect  of  obscurants 
on  rate  of  acquisition  and  rate  of  engagement.  In  Section  II  we  introduce  and  justify  the 
assumptions  behind  Rotman's  smoke  model,  and  show  how  to  extract  the  model 
parameters  for  a  given  obscurant  condition.  In  Section  m  we  write  down  a  formal  solution 
to  the  model,  then  show  that  the  problem  maps  onto  a  well-known  problem  in  the  statistical 
mechanics  of  spin  systems.  The  formalism  which  we  use  to  obtain  solutions  is  somewhat 
nonstandard,  since  the  quantities  that  we  need  to  compute  are  not  the  ones  usually  required 
in  the  physical  problem.  We  proceed  to  obtain  the  analytic  solution  to  our  version  of  the 


^  SJl.  Rotman,  "Modeling  human  search  and  target  acquisition  performance:  n.  Simulating  multiple 
observers  in  dynamic  scenarios,”  Opt.  Eng.  28  (1 1)  (1989). 


problem  in  Section  IV.  In  Section  V  we  pose  and  solve  a  different  problem  within  the 

context  of  Rotman's  model;  namely,  computing  the  probability  of  maintaining  an 

uninterrupted,  unobscured  line  of  sight  for  some  fixed  time  interval.  This  addresses  # 

questions  of  the  success  of  a  subsequent  engagement  phase  in  the  case,  for  example,  of  a 

command  guided  munition.  The  results  of  Sections  IV  and  V  are  summarized  concisely  in 

Section  VI. 
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II.  JUSTIFICATION  OI'  MODEL 


There  are  two  distinct  assumptions  which  underlie  Rotman's  approach  and  which 
must  be  addressed  before  we  consider  the  model  in  detail.  The  Hrst  is  the  two-level 
assumption,  the  second  that  there  is  a  characteristic  time  that  governs  the  transitions 
between  the  two  levels.  In  this  chapter  we  consider  each  of  these  issues  in  turn,  with 
special  attention  to  the  methods  of  calculation  for  the  model  parameters  from  operational 
conditions. 

A.  THE  TWO  LEVELS  OF  OBSCURATION 

The  model  for  smoke  which  Rotman  uses  is  a  two-level  model.  He  assumes  that 
there  is  a  certain  acquisition  rate  that  is  relevant  for  the  "clear"  condition,  and  another 
slower  one  for  the  "snx)ked"  condition.  (In  this  paper  we  specialize  to  the  case  where  the 
slower  rate  is  zero.)  This  is  clearly  a  simplification  of  the  true  situation,  where  any  given 
level  of  obscuration  is  in  principle  possible.  Can  a  model  which  makes  such  an  apparently 
extreme  simplification  have  any  relevance  at  all?  We  contend  that  it  does.  The  reason  for 
this  is  that  target  acquisition  parameters  vary  appreciably  only  over  a  narrow  range  of 
obscuration  conditions.  Most  of  the  time  the  target  visibility  is  either  essentially  un¬ 
degraded  (compared  to  ambient  conditions)  or  essentially  zero.  The  intuitive  justification 
for  this  statement  comes  from  the  commonplace  observation  that  smoke  plumes  and  clouds 
seem  to  have  well  defined  edges,  and  even  holes. 

In  order  to  quantitatively  justify  our  claim,  we  need  to  discuss  the  way  in  which 
smoke  affects  target  visibility.  We  do  so  in  the  context  of  the  Night  Vision  Laboratory's 
Static  ModeP  for  target  acquisition  using  thermal  imaging  sensors.  According  to  this 
noodel,  the  rate  at  which  the  target  is  acquired  is  prcqrortional  to  the  probability  of  ultimately 
detecting  the  target  (Poo),  with  the  constant  of  proportionality  being  0.29  inverse  seconds. 
Poe  is  determined  by  the  ratio  of  the  number  of  resolvable  cycles  on  target  (N)  to  a  cycle 
criterion  constant  which  is  determined  by  the  cognitive  level  of  the  acquisition  (N50)  by  a 


2  J.A.  Rarches,  WJl.  Lawson.  LP.  Obcrt,  RJ.  Bergemann,  T.W.  Cassidy,  and  J.M.  Swenson,  "Night 
Vision  Laboratory  static  performance  model  for  dteimal  viewing  systems,"  Research  and  Devekqxnent 
Tech.  I'-spL  ECOM-7043.  U.  S.  Army  Electronics  Command,  FL  Monmouth.  NJ  (April  1975). 
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target  transfer  probability  function  (TTPF).  N  in  turn  is  computed,  in  units  of  the  angular 
target  size  s,  from  the  apparent  thermal  contrast  of  the  target  (ATa)  by  applying  the  inverse 
of  the  minimum  resolvable  temperature  (MRT)  function,  which  is  the  function  that 
quantitatively  characterizes  the  thermal  sensor.  Finally  ATa  may  be  computed  from  the 
actual  thermal  contrast  (AT)  given  the  attenuation  properties  of  the  obscurant.  We  assume 
that  Beer’s  Law  is  valid;  that  is,  that  the  attenuation  is  exponential  in  the  concentration 
length  ((X)  with  extinction  coefficient  cl 

Before  translating  these  words  into  equations  and  "turning  the  crank,"  it  is 
worthwhile  to  observe  how  it  is  that  the  continuous  variable  CL  can  be  transformed  into  an 
almost  dichotomous  target  acquisition  probability.  There  are  two  links  in  the  foregoing 
chain  of  relationships  that  can  cause  this  to  happen.  First,  Beer's  Law  contains  the  CL  in 
the  exponential,  thereby  amplifying  the  effect  of  small  changes  in  that  variable.  Second, 
the  TTPF  is  a  threshold  function — a  "soft"  one  to  be  sure,  but  when  combined  with  the 
amplification  from  Beer's  Law  the  effect  can  be  essentially  on-off.  Presently  we  will 
derive  the  conditions  under  which  this  is  so. 

Summarizing  these  statements  symbolically,  we  have 

Poo  =  TTPF(5  MRT-1(AT  e-«  CL)/N5o)  ,  (IM) 

where  MRT"1  is  the  functional  inverse  of  the  MRT.  The  TTPF  is  given  in  the  Model  to  be 

^.7+0.7x 

In  order  to  get  an  analytic  result,  we  invoke  a  widely  used  exponential  approximation  to  the 
MRT  function^.^  so  that  it  can  be  inverted  explicitly: 

MRT(v)  =  MRTq  cxp(P^^v)  =»  MRT-kAT)  =  .  01-3) 

where  typical  values  of  the  constants  MRTq  and  are  4  x  10“^  ®C  and  6  x  10“^ 
radians,  respectively.  Note  that  in  this  approximation  the  target  contrast  must  exceed 
MRTo  in  order  to  be  detectable. 


^  R.E.  Roberts,  "A  Simplified  Approach  to  Analyses  of  Infrared  Sensor  Performance  versus  Weather 
Theory  and  Application  to  the  Hanford  Data  Base  (U),"  IDA  Paper  P-1284,  October  1977 
(CONFIDENTIAL). 

^  R.E.  Roberts  and  L.M.  Biberman,  "Impact  of  Atmospheric  Propagation  on  Present  and  Future- 

Generation  Infrared  Imaging  Systems:  A  Reappraisal  of  3-S  pm  FLIR  Performance  (U),’  IDA  P^)^  # 

P-1284.  October  1977  (CONFIDENTIAL). 
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Given  these  relations  we  can 
compute  the  acquisition  probability  as  a 
function  of  CL  for  any  sensor  and  scenario 
conditions  we  choose.  One  set  of  examples 
is  shown  in  Figure  II- 1.  We  have  chosen 
the  sensor  parameters  to  be  as  given  above 
for  all  cases.  The  scenario  dependent 
parameters  for  the  Base  Case  are  AT  =  1®, 
d  =  2  milliradians,  and  N50  =  3.5  cycles. 
Each  Excursion  differs  from  the  Base  Case 
in  one  parameter,  respectively:  (1)  AT  = 
0.5°,  (2)  d  =  1  milliradians,  and  (3)  N50  = 
1  cycle.  Our  intuitive  expectation  is 
confirmed  for  these  cases.  We  observe  the 
transition  from  low  to  high  acquisition 
probability  occurring  over  a  narrow  range 
of  concentradon  length. 


Figure  11-1.  The  value  of  Poo  goes  from 
high  to  low  values  over  a  narrow  range 
of  CL.  The  Base  Case  and  the  various 
excursions  are  described  in  the  text. 


In  order  to  elucidate  the  specific  conditions  where  this  situation  will  occur,  some 
approximations  are  needed.  We  estimate  the  range  of  CL  over  which  Poo  varies 
appreciably  to  be  the  inverse  of  the  derivative  of  Poo  evaluated  at  the  Poo  =  50%  point. 
(Note  that  this  implies  that,  in  the  absence  of  obscurant,  the  target  can  be  acquired  with  at 
least  50%  probability.)  Taking  the  derivative  of  Equation  II- 1  and  using  11-2  and  II-3 
yields 


ACL 


j  IgPsj^ 

ou 
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The  factor  1.18  arises  in  the  derivative  of  the  TTPF.  This  relation  is  less  simple  than  it 
looks.  This  is  because  it  is  not  really  independent  of  CL,  since  we  have  constrained  CL  to 
its  50%  value.  It  is  also  desirable  to  remove  the  "engineering”  parameters  from  the  right 
hand  side  to  the  extent  feasible.  After  some  straightforward  algebra  we  obtain  from 
Equation  11-4  the  condition 


ACL  ,  ,„ln(AT5o/MRTo) 

7s — »  1.18 - 

ln(AT/AT5o) 
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where  CL50  is  the  CL  at  which  the  AT  target  is  50%  detectable,  and  ATsq  is  the  target 
contrast  at  which  it  would  be  50%  detectable  in  the  absence  of  obscurants;  both  definitions 
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can  be  resolved  using  Equations  11- 1  to  11-3.  Now,  the  target  size,  sensor  resolution,  cycle 
criterion  and  smoke  extinction  coefHcient  all  have  been  absorbed  into  these  two  new 
parameters.  Thus  our  mild  approximation 
has  gained  us  some  generality. 


Note  that,  by  hypothesis,  AT  > 
AT50  >  MRTq.  so  both  of  the  logarithms 
are  positive.  The  ratio  of  logarithms  is 
smaller  than  (xie  subject  to  the  condition 


AT> 


(AT5o)2 

MRTo  ’ 
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which  is  the  regime  in  which  our  two-level 
approximation  is  valid;  note  however  that 
due  to  the  logarithmic  figure  of  merit  the 
two-level  hypothesis  degrades  gracefully. 


0.5 

0.4 

0.3p 

0.23 

0.1 
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Given  that  this  cmdition  holds,  one 
needs  to  be  able  to  extract  the  necessary 
model  parameters  from  the  smoke  data. 
These  are  the  value  of  P«»  for  the 
unobscured  condition.  Phi,  and  the  fraction 


Figure  11*2.  Illustration  of  computation 
of  model  parameters  from  plume  CL 
probability  density  function  (PDF) 
and  the  P«,  function. 


of  time  that  is  spent  in  the  unobscured  condition,  q. 


Our  approach  to  determining  i)  is  illustrated  in  Figure  11-2.  We  assume  a  known 
approximate  distribution  function,  f(CL),  which  describes  the  frequency  of  occurrence  of  a 
given  concentration  length  of  obscurant  At  one  level  of  approximation,  we  could  estimate 
T\  by  integrating  the  CL  distribution  up  to  the  point  of  inflection  of  the  Poe  curve.  To 
account  for  the  fact  that  our  "knife"  is  a  somewhat  dull  one,  we  improve  our  precision  by 
con:q)uting  the  integral 

r\  =  JPoo(CL)  f(CL)  dCL  .  (H-?) 


The  factor  in  front  of  the  integral  corrects  for  the  fact  that  Poe  may  not  saturate  at  unity  as 
CL  goes  to  zero.  Similarly  the  appropriate  estimate  for  Phi  is  given  by 


f(CL)  dCL  . 
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A  corresponding  Plow  may  be  computed  (for  comparison  to  zero)  by  replacing  one  power 
of  Poo((X)  in  Equation  II-8  by  the  complementary  quantity  [Poo(O)  -  Poo(CL)]  . 

B  .  THE  CHARACTERISTIC  TIME  SCALE 

A  validated  thecxetical  foundation  for  the  determination  of  a  characteristic  time  scale 
for  the  holes  and  gaps  in  smoke  plumes  has  recently  become  available.  Hoock  and 
Sutherland^  have  applied  the  standard  results  of  Kolmogorov’s  turbulence  theory  to  the 
optical  transmission  problem.  They  begin  by  computing  the  spatial  structure  function  of 
the  CL  through  a  smoke  plume.  The  spatial  correlations  are  then  linked  to  time  by  invoking 
a  "frozen  plume"  approximation.  This  amounts  to  simply  dividing  the  spatial  scale  by  the 
transverse  component  of  the  ambient  wind  velocity. 

The  spatial  structure  function 
presented  in  Hoock  and  Sutherland  is 
shown  as  a  solid  line  in  Figure  11-3.  The 
horizontal  axis  is  length  expressed  in  units 
of  the  plume  thickness,  L  It  is  interesting 
that  the  function  exhibits  a  transition  for 
lengths  equal  to  the  thickness  of  the  plume. 

Fot  shorter  lengths,  the  logarithmic  slope 
approaches  5/3;  for  larger  lengths  it 
approaches  2/3. 

This  curve  actually  represents  a 
special  case;  namely  that  the  "outer  length 
scale,"  Lo,  approaches  infinity.  This  Lq  is 
the  quantity  that  usually  arises  in  turbulence 
problems.  Roughly  speaking,  it  represents 
the  length  scale  at  which  energy  is  injected  into  the  system.  Its  value  typically  is  below  or 
around  100  meters.  To  accommodate  finite  Lot  we  need  to  go  back  to  the  defining  equation 
(Eq.  6.89  on  p.  420  of  Ref.  4)  and  numerically  reevaluate.  The  results  for  particular  values 
of  Lo  are  shown  as  dashed  and  dotted  curves  in  the  figure.  The  net  effect  is  that  the  curve 


Figure  ll-S.  The  smoke  plume  spatial 
structure  function  for  various 
values  of  L. 


^  Donald  W.  Hoock  and  Robert  A.  Sutherland,  "Obscuration  Countermeasures,"  in  The  Irtfrared  and 
Electro-Optical  Handbook,  Vol.  7,  Countenneasure  Systems  (David  Ptdiock,  ed.)  1993. 
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is  unchanged  for  small  values  of  separation,  but  saturates  quickly  as  the  separation 
approaches  the  outer  length  scale. 


As  we  shall  see  later,  the  temporal 
correlations  in  our  model  are  assumed  to  be 
exponential.  So  the  parameter  we  require  is 
the  time  at  which  the  structure  function 
reaches  (1-1/e)  =  0.63  of  its  maximum 
value.  Our  time  scale  is  therefore  of  the 
form 


where  c  is  a  dimensionless  number  of  order 
unity  that  varies  slowly  with  Lo/^.  See 
Figure  11-4.  Given  that  the  uncertainties  in 
determining  Lq  and  Vc  will  certainly  exceed 
10%,  there  is  no  penalty  in  assigning  c  a 
constant  value  of  O.SS. 


Figure  IM.  The  L  dependence  of  the 
"constant”  factor  c. 


( 


III.  STATISTICAL  MECHANICAL  METHOD 

Motivated  by  the  discussion  in  Section  n,  the  probability  of  an  observer  who  will 
eventually  acquire  the  target  of  not  seeing  a  target  by  a  time  t  is  assumed  to  be  governed  by 
a  differential  equation  of  the  form 

^  .  (Ui-i) 

where  a  is  the  usual  target  acquisidon  time  constant  and  Ti(t)  is  a  random  function  of  time 
describing  the  ability  to  see  through  the  obscuration.  When  the  line  of  sight  is  completely 
blocked  il  =  0;  when  the  line  of  sight  is  clear,  =  1.  Intermediate  values  model  partial 
seeing:  this  paper  deals  primarily  with  the  simpler  binary  case  of  complete  visibility  versus 
complete  blockage;  generalizations  to  partial  seeing  will  be  indicated  where  appropriate. 
Eq.  (IIM)  embodies  the  assumption  that  the  important  process  in  target  acquisition  is  the 
total  time  available  with  a  clear  line  of  sight.  This  is  consistent  with  the  modeling  of  the 
acquisition  process  as  a  series  of  short  glimpses  at  different  parts  of  the  field  of  view.^ 

The  solution  to  Eq.  (IIM)  is  immediate. 

t 

-  a  J  Ti(x)  dx 

Pf^l(t)  =  <C  ®  >  .  (ni-2) 

where  the  brackets  (<>)  denote  an  average  over  the  random  obscuration  T|(t).  In  the 
unobscured  limit  for  which  11(0  =  1  for  all  values  of  t,  Pfeii  (t)  =  xp(-  at)  and  therefore  the 
probability  of  successful  target  acquisition  for  this  observer  is  Psijccecd(0  =  1  -  cxp(-  at). 

This  expectation  value  of  an  exponential  has  a  strong  connection  to  statistical 
models  of  ferromagnetism.  To  develop  this  analogy,  a  "spin"  variable  s(t)  is  defined  by 
11(0  =  [1  +  s(t)]/2  so  that  s(t)  =  +1  and  -1  correspond  to  il(t)  =  1,  and  0,  respectively. 
Thus,  the  expectation  value  in  Eq.  (111-2)  is  rewritten  as 


An  alternative  hypothesis  would  be  to  require  a  fixed  amount  of  continuous  uninterrupted  observation. 


This  will  be  addressed  in  Section  V. 
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-  a  J  ii(x) 


>=  t°^<C 


-  (o/2)Js(x)dx 


(ni-3) 


It  is  the  average  over  the  spin  variable  that  will  be  calculated  explicitly.  In  order  to  apply 
the  standard  magnedc  models,  it  is  convenient  to  break  up  the  integral  over  time  into  a  sum. 
Divide  the  interval  from  0  to  t  into  N  segments  of  duration  A,  t  =  NA;  and  rewrite  the 
integral  as  a  sum. 


f  N 

J  s(x)dx  =  A^  Sj 


(ni-4) 


The  large  N  (equivalendy,  small  A)  limit  will  be  taken  at  the  end  of  the  calculation. 


The  problem  has  now  been  recast  into  the  fonn  of  a  one-dimensional  chain  of  spins 
taking  on  the  values  of  ±  1  (referred  to  as  "spin  up"  and  "spin  down"  states);  such  a  model 
is  termed  an  Ising  model.  Denoted  by  P({s})  the  probability  that  the  set  of  N+1  spins  (s) 
takes  on  a  pardcular  arrangement  of  up  and  down  spins,  the  average  over  spins  required  for 
Eq.  (ni-3)  can  be  written  as: 


.(oA/2)i  S:  -  (“A/2)Ssi 

J  XP((s))c 

(s) 


(ni-5) 


This  can  be  placed  into  more  usual  thermodynamic  form  by  representing  P({s})  in  terms  of 
a  Hamiltonian  or  energy  function^ 


P({S))  occ 


H((s))  +  hXsi 


(ni-6) 


where  the  term  linear  in  the  sum  of  the  spin  values  has  been  separated  from  the  remainder 
of  the  Hamiltonian,  H({s}).  The  coefficient  of  the  linear  term,  h,  represents  the  magnetic 
field  of  the  spin  system  and  serves  to  control  the  average  value  of  the  spin.  Thus  a  large 
positive  h  will  induce  the  spins  to  point  up;  translating  back  to  the  obscurant  model,  a  large 
positive  h  increases  the  chance  of  a  clear  line  of  sight. 


To  be  precise,  in  a  true  thermodynaniics  system  this  dimensionless  Hamiltonian  would  correspond  to 
the  usual  Hamiltonian  divided  by  kT  where  k  is  Boltzmann’s  constant  and  T  the  absolute  temperature. 
The  same  is  tnie  of  the  dimensionless  Gibbs  hee  energy  given  in  Eq.  Gn-8). 


Inserting  this  expression  for  P({s} }  into  Eq.  (III-5),  the  first  result  of  this  section  is 
obtained. 


.(aA/2) 

<e 


N  n 

H((s))  +  hSsi-(oA/2)X  s. 
M) 


N 
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S' 


-  H({S))  -hh  £  s. 


i-0 
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Note  that  the  numerator  and  dentxninator  of  this  fraction  is  of  the  same  form  and  differ  only 
in  the  effective  value  of  the  "magnetic  field,"  h.  In  fact,  the  denominator  is  just  the  partition 
function  for  this  system,  Z(h),  and  the  numerator  is  the  same  partition  function  with  a 
modified  argument,  Z(h  -  aA/2).  Using  the  fact  that  the  Gibbs  free  energy,  G,  of  a  system 


is  related  to  the  partition  function  by  Z  =  exp(-  G),  Eq.  (111-7)  can  be  rewritten  as 
^^■(ttA/2)^s.^  Z(h-aA/2)  _  ^G(h)  -  G(h  -  aA/2)  . 


(in-8) 


The  original  problem  of  the  expected  value  of  the  probability  of  not  seeing  a  target  in  the 
presence  of  obscurants  is  thus  shown  to  be  equivalent  to  finding  the  Gibbs  ^e  energy  of  a 
related  Ising  magnetic  system.  This  "smoke  scholium"  allows  the  application  of  a  wide 
range  of  one-dimension  studies  of  Ising  systems  that  have  been  solved  for  a  number  of 
different  Hamiltonians.  Each  choice  of  Hamiltonian  corresponds  to  a  different  distribution 
of  smoke  properties;  a  specific  choice  will  have  to  be  made  to  produce  a  complete 
descripticm  of  the  acquisition  problem. 


Some  results  can  be  obtained  directly  from  the  thermodynamic  formalism  and  can 
be  used  to  check  the  results  of  any  detailed  calculation. 

(1)  For  a  «  1,  G(h)-G(h-(xA/2) »  aA/2N<s>.  This  is  equivalent  to; 


I 

-  a  J  -n(x)dx 

<e  “  >  -  .  (ni-9a) 

That  is,  for  small  a,  the  time  constant  is  reduced  by  the  average  probability  of  obscuration. 

(2)  For  t «  Tc,  where  Tc  is  the  smdce  correlation  time. 


-  a  J  Ti(x)dx 


<e 


>  »  p^e”®*  +  p. 


(ni-9b) 
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where  p+  and  p.  are  the  probabilities  that  the  initial  state  of  the  system  was  initially 
(at  t  =  0)  in  the  spin  up  (clear  vision)  state  and  spin  down  (obscured  vision)  state, 
respectively. 

(3)  For  t  or  N  very  large,  the  Gibbs  free  energy  should  be  "intensive";  that  is,  G 
should  be  proportional  to  N.  In  that  case, 

t 

-  a  J 

<e  “  >  »  Ce‘“‘"‘  .  (ni-9c) 

where  C  and  Oeff  are  constants. 

In  Section  IV,  a  particular  form  of  the  Hamiltonian  will  be  chosen  that  allows  for  a 
simple  calculation  of  the  complete  result. 
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IV.  SMOKE  MODEL  RESULTS 


The  complete  solution  to  Eq.  (111-2)  can  be  obtained  easily  if  a  simple  form  of  the 
Hamiltonian  is  chosen.  The  one  chosen  is  the  simplest  Hamiltonian  to  solve  and  also 
corresponds  (as  will  be  shown)  to  the  assumption  of  the  Kowalczyk-Rotman  Monte  Carlo 
approach.  The  nearest-neighbor  Ising  model  contains  interactions  only  between  adjacent 
spins: 

H=-e£(SjS-ll  (IV-1) 

i=j±i 


where  e  is  an  "interaction  energy"  and  serves  to  control  the  coupling  or  correlation  between 
spins.  A  large  value  of  this  interaction  will  produce  a  strong  correlation  between  spins 
corresponding  to  longer  obscurant  correlation  tin^s.  In  this  section  the  partition  function 
ratio  for  the  Ising  model  Hamiltonian  will  first  be  calculated  using  conventional  statistical 
mechanical  techniques.  Subsection  A  will  review  the  solution  methods  for  the  usual  Ising 
model;  subsection  B  will  derive  the  results  for  the  partition  function  ratio  which  describes 
the  probability  of  not  seeing,  Eq.  (111-8).  The  result  shows  that  the  choice  of  the  Ising 
Hamiltonian  is  equivalent  to  assuming  exponential  distributions  of  both  the  clear  and  block 
intervals.  This  being  the  case,  a  Markov  random  field  approach  is  appropriate.  The 
rederivation  of  the  results  using  the  Maikov  aj^iroach  will  be  given  in  subsection  C. 


A.  ISING  MODEL  SOLUTION 

The  Ising  model  Hamiltonian,  Eq.  (IV-l),  is  solved  conventionally  by  introducing 
the  transfer  matrix,  Tij  =  exp  (-Hij)  where  i  and  j  indicate  the  states  (+1  or  -1)  of  adjacent 
spins.  Including  the  magnetic  field  contribution,  T  is  given  by 


T  = 


.-2£ 


e-2«  e  -** 


(IV-2) 


Now  imagine  M  Ising  spins  arranged  not  in  a  straight  chain  but  in  a  circle  so  that  the  last 
spin  and  first  spin  are  identical  (si  =  sm)-  Then  the  partition  function  which  is  the  sum 
over  all  states  can  be  represented  by 

Z  =  traceT^  .  (IV-3) 
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The  matrix  product  takes  care  of  the  sum  over  the  states  of  the  spins  and  the  trace  is  the 
representation  of  the  periodic  boundary  condition  imposed.*  The  trace  can  be  evaluated  by 
recognizing  that  T  can  be  diagonalized  with  some  orthogonal  matrix,  T  =  ODO“^  where 


D  = 


is  the  matrix  of  eigenvalues  of  T  (X+  >  L): 


0 


0 

X. 


av-4) 


X^=  cosh(h)  ±  J  sinh^(h)  +  e“ 


(IV-5a) 


and 


O  = 


COS0  sin6 
-sin0  COS0 


(IV-5b) 


where  0  is  the  rotation  angle  of  the  orthogonal  matrix.  Then  the  partition  function  is 


Z= 


*•  A,^ ,  as  M 


av-6) 


A  number  of  quantities  can  be  calculated  in  addition  to  the  partition  function.  The 

average  value  of  the  spin  <s>  is  given  by 

2/a  •  sinh(h) 

<s>  =  cos  0  -  sin  0  = 


sinh  (h)  +  e 


-4c 


av-7) 


The  probabilities  of  a  particular  site  being  spin  up  (clear  line  of  sight)  or  spin  down 
(blocked  line  of  sight)  were  denoted  as  p.*.  and  p-  in  Sec.  m  and  are  given  by  p^.  =  cos^ 
and  p-  =  sin^,  respectively.  Hnally,  the  spin-^in  correlation  function  is  given  by 


<Sj  Sj.  >  =  <s>^  +  (1  -  <s>^)  (^ )  *  ^  * 


(IV-8) 


*  The  matrix  itself  is  a  matrix  of  partition  functions  for  fixed  boundary  conditions  of  the  open  chain 

of  M-f  1  spins,  the  -h-  element  of  which  corresponds  to  the  partition  function  with  sq  and  Sm  both  iq),  # 

etc. 
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This  completes  the  solution  of  the  usual  Ising  model  for  the  infmite  chain,^  but  is  not 
sufficient  evaluate  the  ratio  of  partition  functions  in  Eq.  (111-8)  that  require  consideration 
of  finite  chain  effects.  For  future  use.  the  appropriate  limits  of  the  infinite  chain  Ising 
quantities  for  the  small  A  limit  are  needed.  Consider  a  scaled  magnetic  field  h  =  /i  A,  and 
near-neighbor  interaction  parameter,  exp  (-2e)  =  £  A.  Then 

=  1  ±  A|3/2 

^ 


<s> 


2h 

P 


<s(o)s(t)>  =  <s>^  +  (1  -  <s>^)  . 


av-9) 


These  changes  of  scale  were  chosen  to  hold  the  average  value  of  the  spin,  <s>  (related  to 
the  average  probability  of  obscurant),  and  the  correlation  time  of  the  obscurant  (given  by 
1/^)  fixed  in  the  limit  A  — >  0. 


B .  APPLICATION  TO  THE  SMOKE  SCHOLIUM 

Instead  of  a  circle  of  M  spins  with  a  single  interaction  between  all  the  spin  pairs, 
imagine  that  M-N  spins  interact  with  the  Ising  Hamiltonian  and  magnetic  Held  h  with 
transfer  matrix  T,  while  the  remaining  spins  have  a  field  h-aA/2  and  transfer  matrix,  Tq. 
Then  the  partition  function  for  this  case,  denoted  Z(a),  is  simply 

Z(a)  =  trace  (TM-N  To^  )  .  (TV-IO) 

The  desired  ratio  of  partition  functions  in  Eq.  (111-8)  is  Z(a)/Z(0).  Denoting  by  D  and  Dq 
the  eigenvalue  matrices  for  the  two  cases,  and  O  and  Oa  the  corresponding  rotation 
matrices,  this  is: 

Z(a)  *  trace  (DM-No-lOoP„  No„-lO)  .  (IV-1 1) 

The  product  of  the  two  rotation  matrices  is  itself  a  rotation  matrix  with  angle  56  =  0-6a 
Taking  the  large  M  and  large  N  limits  (with  NA  ^  t  Hxed),  the  partition  function  can  be 
straightforwardly  evaluated.  Taking  the  ratio  of  Z(a)/Z(0)  and  replacing  the  factor  of 
exp(-oet),  the  final  expression  for  Eq.  (111-2)  is 


^  It  is  easy  to  show  that  the  use  of  periodic  boundary  conditions  is  not  essential  in  the  solution. 


IV-3 


-  o  I  Tl(x)dx 

<e  “  >=  e'°  '^*''’"(cos^8e  +  sin26ee‘'^‘]  .  aV-12a) 


where 


Pa 


cos^86 


=  (3  -  a<s>  (1  -  <s>^) 

P  -  a  <s> 


1  + 


Pa 


av-12b) 


sin^50 


,  P  -  a  <s> 
Pa 


The  result  is  a  mixture  of  two  exponentials  whose  time  constants  are  functions  of  the  target 
acquisition  time  constant,  orl,  the  smoke  ctxrelation  time  P~^  and  the  average  probability 
of  clear  line  of  sight,  <ti>.  <n>  is  given  by  (1  +  <s>)/2.  Rgure  FV-l  shows  a  comparison 
of  the  line  of  sight  probabilities  for  a  number  of  cases.  The  figure  shows  the  unobscured 
result  (<s>  =  1)  and  two  variants  of  the  obscured  cases  for  <s>  =  0  and  <s>  =  -  0.5 
(corresponding  to  <cn>  =  1,  <n>  =  0.5,  and  <n>  =  0.25,  respectively)  and  a  range  of 
values  of  p/a;  p/a  =  0.5, 1.0, 2.0,  and  5.0.  Recall  that  <n>  is  the  average  probability  of  a 
clear  line  of  sight  while  P  measures  the  scale  of  the  clouds  and  gaps  (large  p  being 
equivalent  to  short  intervals  of  clouds  and  short  gaps  for  fixed  <ti>). 

Note  that  the  curves  for  each  value  of  <ii>  are  bunched  together  even  though  the 
values  of  p/a  vary  by  a  factor  of  25.  For  a  fixed  value  of  <t|>,  the  probabilities  are 
bounded  by  the  P  =  0  and  p  =  <»  limits; 

^  ”  oiXl-e") 

which  agrees  with  Eqs.  (III-9a>-(ni-9b)  (identifying  p+  =  and  p.  =  1  The 

relationship  between  the  exact  result  and  the  bounds^^  is  illustrated  in  Fig  rV-2,  which 
reproduces  the  <T|>  =  0.5  cases  shown  in  Fig.  FV-l  together  with  the  bounds. 


10 


In  fact,  the  P  =  0  limit  is  a  craitct  approximation  to  second  order  in  t;  P  dependence  only  enters  for  t^  # 

and  higher  terms. 


IV-4 


Probability 


Figure  IV>1.  Search  probabilities  as  a  function  of  a  t  for 
unobscured  and  obscured  cases. 


IV-5 


Probability 


Search  Timelines 


P^a*5.0:  <7)>sOJ 
p/a  s  2.0:  <ii>s0.5 
P/Osl.0;  <Tp>sO.S 
P/asOJ;  <11>a05 

P/o  =  0.0:  crj>  =  03 


Tima 


•iw  ; '  •»'  «"obtcuiM 


Hnally,  fw  t  large,  the  last  limit  described  in  (in-9)  can  be  recovered: 


cos^se 


(IV- 14b) 


so  that  Oeff  =  (a  -  ^  +  P)/2. 

The  final  result,  Eq.  (rV-12),  reproduces  analytically  the  Monte  Carlo  results  of 
Kowalczyk-Rotman.  The  essential  similarity  connection  is  that  Kowalczyk-Rotman 
assume  that  the  smoke  obscurations  and  the  spaces  between  them  both  have  expcmential 
distributions.  In  this  Ising  approach,  one  requires  the  probability  that  if  a  spin  at  some 
time,  say  t  =:  0,  is  up,  that  all  the  spins  are  up  for  times  between  0  and  t.  This  can  be 
calculated  using  the  transfer  matrix  approach  by  inserting  into  the  trace  of  Eq.  (IV-3)  the 
product  of  n  projection  operators  (n  A  =  t )  of  the  form 

1  0 
"  0  0 


and  then  taking  the  limit  of  small  A.  Similarly,  using  the  complementary  projection 
operator  P.  =  1  -  P.f,  one  obtains  the  probability  of  staying  in  a  spin  down  state  (line  of 
sight  blocked).  The  resulting  probabilities  are  exponentials: 

P±(t)=  e^‘  (IV-15a) 


K 


=  P  •  <s>) 
2 


(IV-15b) 


P.  =  P  . 

The  two  exponents  ^  could  be  used  to  characterize  the  system  instead  of  <n>  and  p. 


C.  MARKOV  RANDOM  HELD  APPROACH 

The  Ising  model  Hamiltonian  is  therefne  equivalent  to  a  two  exponent  description 
of  the  intervals  of  blocked  and  clear  line  of  sight  There  is  an  alternative  derivation  of  the 
result  given  in  Eq.  (TV- 12)  which  exploits  that  relationship  directly.  This  uses,  instead  of 
the  statistical  mechanical  analogy,  the  language  of  Markov  processes.^^  To  apply  the 


This  method  was  suggested  to  the  authors  by  Dr.  AmnonDalcher  of  IDA 


Markov  method,  consider  a  three-state  system;  (1)  the  spin  down  state  corresponding  to 
blocked  sight;  (2)  the  spin  up  state  (clear  line  of  sight);  and  (3)  a  "done"  state, 
corresponding  to  having  acquired  the  target.  A  transition  matrix  between  these  states  can 
be  written  down  immediately  in  tenns  of  a  and  P±.  Normally,  a  three-state  system  would 
require  a  3  x  3  transition  matrix  for  a  complete  description;  however,  in  this  case,  the  done 
state  has  no  transitions  to  either  the  up  or  down  states  and  therefore  a  transfer  matrix  to 
describe  just  the  transition  of  the  up  and  down  states  suffices. 


T 


-(P^  +  a) 

P.  -P- 


(IV-16) 


Thus,  the  up  state  decays  through  two  channels;  (1)  at  a  rate  p.^.  into  the  down  channel  and 
(2)  at  a  rate  of  a  into  the  done  channel.  The  down  state  can  only  transition  into  the  up 
channel  and  does  so  at  a  rate  of  P_  The  exponential  of  T  gives  the  transition  probabilities; 


-  n  -  c^  ‘(r  -X*)] 

Pa 

where  X*  are  the  eigenvalues  used  in  Eq.  (IV- 12); 

,±  (P  +  a)  .  Pa 

X  =  — y-±Y  . 


(IV-17) 


(IV- 18) 


This  can  be  used  to  describe  the  probability  of  being  in  the  up  and  down  states  for  any 
initial  distribution  of  states.  The  probability  of  not  seeing  the  target  is  just  the  probability  of 
being  in  either  the  up  or  down  state  (since  otherwise  the  system  is  in  the  done  state).  For 
any  initial  distribution  of  up  and  down  states  (v+,  v_),  with  v+  +  v_  =  1,  the  probability  of 
not  seeing  is  therefore^^; 


‘  ( -  X’  -  v^  a)  -  e^'^  ( -  X^  -  v^a) 


(IV-19) 


For  the  a  priori  estimate  of  the  distribution  of  states,  v+  =  p+  =  (1  +  <s>)/2,  this  reduces  to 
Eq.  (TV- 12).  On  the  other  hand,  if  one  knows  the  state  at  t  =  0  to  be  definitely  clear  or 
blocked  the  appropriate  choice  of  v+  is  1  ot  0  respectively. 


^  2  This  means  that  the  transition  matrix  given  here  does  not  "conserve  probability."  The  lost  {Hobability 
corre^nds  to  populating  the  done  state. 

^  ^  This  result  is  more  general  than  that  given  in  Eq.  (TV- 12)  using  the  statistical  mechanical  formulation; 

however,  the  same  general  result  can  be  obtained  in  the  partition  function  approach  by  using  fixed  # 

rather  than  free  boundary  conditions. 
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Generalizations  of  either  the  partition  function  approach  or  the  Markov  approach 
could  involve  the  introduction  of  intermediate  states  between  clear  and  totally  blocked  lines 
of  sight.  In  addition,  more  complex  interactions  can  be  modeled  in  the  Hamiltonian  to 
provide  something  other  than  simple  exponential  distributions  of  the  clear  and  blocked 
intervals.  However,  any  of  these  generalizations  involves  the  introduction  of  additional 
parameters.  The  experimental  data  available  at  this  time  does  not  indicate  the  need  for  any 
such  extension. 


IV-9 


V.  MODELING  EXTENDED  INTERVALS 


A.  APPLICATION  OF  THE  ISING-MARKOV  MODELS  TO 
TIME  INTERVALS 

In  the  preceding  sections,  the  target  acquisition  task  was  modeled  as  a  cumulative 
process  for  which  the  total  time  available  on  target  was  the  relevant  quantity,  independent 
of  how  that  time  was  accumulated.  This  is  consonant  with  the  assumption  that  target 
acquisition  is  the  consequence  of  a  series  of  independent  glimpses,  each  of  which  is  of 
short  duration  compared  to  any  change  in  the  visibility  conditions.  Thus,  the  task  is 
performed  equally  well  when  10  seconds  of  clear  viewing  is  provided  whether  it  is 
provided  in  a  single  interval  of  10  second  or  in  10  intervals  of  1  second  each.  Another 
point  of  view  would  be  to  assume  that  a  single  continuous  block  of  time  is  necessary  to 
perform  a  task.  This  would  be  the  case,  for  example,  for  a  connmand-guided  missile,  but 
may  also  be  relevant  for  other  aspects  of  the  target  recognition  process. 

One  can  describe  this  problem  using  the  Ising-Markov  models  of  the  obscurant 
developed  in  Sections  HI  and  IV.  The  line  of  sight  is  alternately  clear  and  blocked.  It  is 
assumed  that  the  task  is  completed  when  a  clear  line  of  sight  is  made  available  for  a  time 
greater  than  or  equal  to  the  required  time  for  the  task,  x.  As  before,  it  is  easy  to  calculate 
the  probability  of  failing,  that  is,  the  probability  as  a  function  of  time,  t,  that  no  such 
interval  of  length  x  or  greater  exists. 

Assume  for  the  moment  that  at  t  =  0,  the  system  is  known  to  be  in  the  spin  up 
(clear)  state.  (Consider  a  series  of  flips  from  spin  up  to  spin  down  and  back  again  occurring 
at  intervals  of  to,  ti,  t2, .  The  total  time  elapsed,  t,  is  the  sum  of  these  intervals. 

Stj  =  t  .  (V-1) 

i 

To  each  transition  between  states,  one  can  associate  a  transition  probability  density  of  e 
^tj  where  tj  is  the  length  of  the  interval  and  the  upper  sign  is  used  for  the  transition  firom 
an  up  state  and  the  lower  for  a  transition  from  a  down  state.  Consider,  for  example  a 
sequence  including  M  flips,  starting  initially  in  the  up  (+)  state.  Then  the  probability 
density  that  the  M  flips  representing  intervals  (  }  will  occur  in  the  total  time  t  is  given  by 


P(M,t,+)  = 


/  V  Mmod  2 


P- 


[M/2) 


-B  t  m-[m/21  .p^,  5(yt.-t) 

^  P*'  <i4i  P  P-'  '  “Vi 


where  (M/2]  is  the  integer  part  of  M/2.  The  corresponding  probability  is  given  by 
integrating  over  the  allowed  values  of  the  sizes  of  the  intervals.  Suppose  that  probability 
that  no  (+)  state  lasts  longer  than  time  t+  and  that  no  minus  (-)  state  lasts  longer  that  t_. 
Then  the  integrals  may  be  performed  if  the  delta  function  is  replaced  by  its  Fourier 
transform: 


r  dco 

J27cP^ 

[M/2]-t-l  M  •  [M/2] 

R- 


(V-3a) 


where  the  factors  R±  are  given  by: 

R±  = 


-l-iO) 


(V-3b) 


Summing  over  all  values  of  M,  one  has 


[P.  +  p-r;') 


(V^) 


The  corresponding  formula  for  beginning  in  the  down  state  is  given  by  interchanging  plus 
and  minus  subscripts.  (If  the  initial  state  is  not  known,  then  one  can  sum  over  the  two 
cases  using  the  a  priori  estimates  of  being  in  the  up  and  down  states  p.|.  =  P-  /p, 
p_  =  P+  /p,  where,  as  before,  P  =  P+  +  . ) 

In  the  case  of  interest,  the  amount  of  time  in  the  down  or  (-)  state  is  unimportant 
(t_  =  oo);  writing  T+  =  x  fOT  simplicity,  Eq.  (V-4)  then  gives  the  probability  that  in  a  time  t, 
beginning  in  the  (+)  state,  no  interval  of  length  x  or  longer  exists  (that  is,  the  probability  of 
failure),  as  follows; 

1.^  Qw  1  -(P++io>>rx 

Ptt  +  X  ^  e’®*  (o)  -  tP)(  1  -  e _ ^ 


(0(ta-iP)  -  p^pL-e 


(V-5a) 


For  starting  in  the  down  (-)  state: 


P(t.-,T)=  0'““- 


«.-iP  +  ip.e^-"‘“>' 


(0(«)-iP)-P,p.e-»'"‘">' 


(V-5b) 


V-2 


The  contour  integrals  in  Eq.  (V-5)  run  along  and  just  below  the  real  to  axis. 
Although  these  integrals  provide  the  exact  results  needed,  the  complicated  form  of  the 
denominators  in  these  integrals  makes  them  difficult  to  evaluate.  One  approach  is  to 
expand  the  denominators  in  powers  of  the  exponential  terms  and  integrate  each  term 
separately.  Note  that  a  term  proportional  to  the  nth  power  of  the  exponential  carries  an  (O 

dependence  of  e“^‘ '  ;  for  t  <  nt,  the  contour  can  be  closed  down  in  the  complex  (o  plane 

and  there  is  no  contribution  to  the  integral.  For  t  >  nx,  the  contour  is  closed  up  and  the 
contributions  from  the  poles  at  to  =  0  and  O)  =  iP  can  be  easily  calculated.*^  Straight¬ 
forward,  but  extremely  tedious,  algebra  gives: 

-Mr 


i  r 


(IV-6a) 


Dj(  j ,  PT. )  -  Dj(  j ,  PT^j )  DjC  j .  PX. )  e 

For  starting  in  the  down  (-)  state: 


-PT  -PT  -  P,T 

^-D2(j,pT.^,)e 


P 


(V-6b) 


P- 


-px 


D,(j,PT.)  +  D^(j.PT.)  -  ^ 


where  Tj  is  a  shifted  time  variable: 


T.  =  t-jt  ;  t>jT 


=  0 


;  t<jT  . 


(V-6C) 


It  is  important  to  note  that  terms  containing  Tj  contribute  only  if  Tj  >  0;  that  is,  a  term  such 
as  Di  (j,  pXj)  carries  with  it  an  implicit  step  function  0(Tj).  The  functions  Di  arc  given  by: 


Note  that  these  are  not  poles  of  the  integrand  before  the  expansion  in  powers  of  the  exponential.  A 
little  thought  shows  that  these  expansions  do  make  sense.  For  example,  one  can  expand  the  simple 
#  integrand  l/((a  -  ooo)  in  powers  of  too.  changing  the  simple  pole  at  a>o  to  a  series  of  higher  order 

poles  at  0)  =  0. 


V-3 


(V-6d) 


=  1 


,  j  =0. 

(V-6e) 

.  j^o  . 

(V-6f) 

This  result,  although  intricate  in  form,  is  relatively  simple  to  describe.  One  can 
describe  the  process  of  waiting  for  a  clear  interval  as  follows.  Starting  initially  in  the  clear 

-P.T 

(up  or  +)  state,  the  probability  that  the  line  of  sight  stays  clear  for  a  time  t  is  e  .  If  the 
line  of  sight  did  not  stay  clear  for  the  whole  interval,  it  must  have  flipped  to  the  blocked 
(down  or  -)  state  at  some  time  before  x.  One  then  must  wait  for  the  blocked  state  to  flip 
and  then  begin  the  next  clear  interval.  If  the  system  flips  before  the  time  x,  then  another 
cycle  of  waiting  for  the  system  to  flip  back  to  the  (*«-)  state  and  then  counting  down  to  x 
begins.  This  alternation  between  waiting  fOT  the  system  to  move  from  the  (-)  to  (+)  states 
and  then  testing  the  duration  of  the  (4-)  state  is  responsible  for  the  introduction  of  the  step 
functions. 

Although  the  complete  expressions  given  in  Eq.  (V-6)  are  daunting,  the  resulting 
probabilities  are  not  themselves  complicated.  Figure  V-1  illustrates  the  probability  of 
success  (having  a  clear  line  of  sight  for  an  interval  greater  than  or  equal  to  x)  for  ^^.x  -  0.5 
and  ^x  =  0.5.  That  is,  the  clear  and  obscured  intervals  both  have  an  average  value  of  2x. 
Two  curves  are  shown;  the  higher  curve,  labeled  P+,  corresponds  to  the  probability  of 
success  starting  from  the  +  (unobscured)  state;  the  lower  is  the  corresponding  probability 
starting  in  the  -  (blocked)  state.  Both  P+  and  P-  are  identically  zero  for  t  /x  <  1,  since  one 
cannot  have  an  interval  of  length  x  until  at  least  that  length  of  time  has  passed.  P+  then 

-B  X 

jumps  discontinuously  to  a  value  of  1  -  e  ,(»  0*61)  corresponding  to  the  probability 
that  the  line  of  sight  has  stayed  continuously  clear.  On  the  other  hand,  P-  is  continuous  at 
t  /x  =  1  since  at  least  one  flip  (from  -  to  +)  must  happen  and  therefore  the  probability 
increases  smoothly  from  0  at  t  /c  =  1.  Note  that  by  t  /x  =  4,  the  probability  of  an  interval  of 
length  at  least  x  is  greater  than  0.8.  For  long  times,  the  distinction  between  P+  and  P- 
becomes  less  important. 


Probability  of  Success 


Figure  V*l.  Probability  of  having  a  clear  Interval  of  length  ^  t 
as  a  function  of  t/t.  (P-t  s  0.5  and  s  0.5) 


Figure  V-2  shows  the  same  curves  for  p_t  =  0.5  and  p+t  =  1.0.  For  these  values, 
the  average  duration  of  a  clear  interval  is  t.  There  is  a  drop  in  the  initial  probability  at  t  / 1  = 
1  to  0.31  (in  comparison  to  the  previous  case)  and  the  rate  of  increase  of  probability  is 
slower.  For  0.8  probability  of  success,  one  must  have  t  /t  >  9. 

Figure  V-3  shows  the  results  for  P-X  =  0.5  and  p+x  =  2.0.  In  this  case  clear 
intervals  last  only  x  /2  on  average;  in  consequence,  the  probability  of  success  is  less  than 
80%  for  t/x  <  25. 


Aside  from  the  first  discontinuity  in  P-f  at  t/T  =1,  the  curves  are  continuous.  The 
step  functions  which  appear  in  Eq.  V-6  only  introduce  discontinuities  in  the  higher  order 
derivatives.  Expanding  the  appropriate  terms  near  the  boundaries,  one  has: 

Dj  (j,z)  +  D2(j,2)c*  « 


D3(j,z)-D3(j,-z)e’*  z«l,j>0 


(V-7) 


so  that  only  the  (j+1)  th  derivatives  exhibit  the  discontinuity. 


Figure  V*2.  Probability  of  having  a  clear  Interval  of  length  ^  t 
as  a  function  of  t/t.  (^-t  =  0.5  and  s  1.0) 


Probability  of  Success 


Figure  V>3.  Probability  of  having  a  clear  Interval  of  length  2  x 

as  a  function  of  t/x.  (P-x  s  0.5  and  p+x  s  2.0)  * 
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B .  SINGLE  POLE  APPROXIMATION 


» 


i 


» 


% 


Although  the  details  of  these  results  depend  on  the  precise  values  of  |Lt  and  P+x  , 
it  is  clear  that  the  overall  structure  is  simple,  particularly  for  longer  times.  The  structure  of 
the  contour  integral  leads  one  to  expect  that,  asymptotically,  these  probabilities  will  be 
dominated  by  a  simple  exponential  form.  In  fact,  these  functions  are  well  represented  by  a 
single  pole  approximation  which  looks  for  the  smallest  root  of  the  denominators  in 
Eqs.  (V-5).  Writing  (o  =  i  K,  the  transcendental  equation  that  needs  to  be  solved  is 

K(K-p)  +  p^P-e~^^*"’'^^  =  0  .  (V-8) 


There  are  two  real  roots  for  k,  the  remaining  (infinitely  many)  roots  are  complex.  One  root 
of  this  equation  is  always  k  =  however,  an  examination  of  Eq.  (V-5)  shows  that  both 
integrals  have  zero  residue  at  that  pole  and  do  not  contribute  to  the  result.  The  other  root 
must  be  determined  numerically  and  the  residue  at  that  pole  gives  the  asymptotic  behavior 
desired.  Denoting  that  root  simply  by  ic,  the  single  pole  approximation  is 


P+  =  1 


P-  =  1 


K  (K  -  p) 

(K-P)(l-t-  ^-o---) 

_ _  e-^ 

2k-P-xk(k-P) 

(K-P)(l-  — ) 
_ P-t-  g-ia 

2k-P-xk(k-P) 


(V-9) 


In  the  special  case  of  P+  -  P_  -(-  P+P-X  =  0,  the  two  roots  are  degenerate  and  the  correct 
asymptotic  form  is  given  by 


P+  =  1  - 


P-  =  1  - 


P_x 


1  +P^P-x2/2 


(V-10) 


l+p_x  -P+t 


l  +  P^P-X^/Z 


Table  V-1  gives  the  values  of  the  transcendental  root  for  a  number  of  values  of  P+x  and 
P-x. 
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Table  V-1.  Values  of  the  Dominant  Pole 


0.1 


.2 


.5 


1.0 


2. 


5. 


0.1 


0.0901 


.2 


0.0811 


1 

2 

5 

0.0358 

0.0131 

0.007 

0.698 

0.0255 

0.0013 

0.1615 

0.0588 

0.0031 

0.2857 

0.1037 

0.0057 

0.4568 

0.1669 

0.0098 

0.6935 

0.2606 

0.0172 

0.1792  0.1608  I  0.1169  0.698 


0.4404 


0.8590  0.7355 


1.5499  1.2798  0.8209  0.4568 


2.6634  2.0683  1.2586  0.6935 


The  single  pole  approximation  is  actually  very  good  for  even  moderate  values  of 
t  /  X  .  Table  V-2  gives  the  errors  for  t  /  x  =  0,  1,  and  2  for  the  cases  given  in  Figs.  V-1, 
V-2,  and  V-3.  The  errors  are  large  at  t/x  =  0,  moderate  at  t/x  =  1  and  negligible  at 
t  /  X  =  2.  This  indicates  that  the  next  (complex)  root  of  Eq.  (V-8)  is  usually  large.  An 
estimate  of  the  location  of  that  root  has  not  been  determined  analytically  but  numerical 
searches  show  it  to  be  typically  10  times  greater  than  k.  Numerical  experiments  have  not 
found  any  nontrivial  cases  (i.e.,  any  case  for  which  the  probabilities  are  not  essentially  0  or 
1)  for  which  the  single  pole  approximation  is  not  essentially  perfect  for  t  /  x  >  2. 

Table  V*2.  Comparison  of  Exact  Result  and  Single  Pole  Approximation 


Because  the  single  pole  approximation  is  so  good,  a  practical  use  of  the  exact  result 
would  be  to  employ  the  exact  result  for  small  t  / 1  (less  than  2  or  3)  and  then  use  the  single 
pole  approximation.  For  convenience  for  such  applicaticms,  the  first  few  Dk(jt  z)  are  given 
below: 

Dj(0,z)  =  1 


D,  (l.z)  =  1  -  z 
Dj(2.z)  =  i[6  -  4z  +  z^] 

(V-lla) 

Dj(3.z)  =  |[60  -  36z  +  9z2 

( 0 .  z)  =0 

D2(1.z)  =  -1 


0^(2. z)  =  .(3  +  z) 

D^O.z)  =  -1(20 +8z  +  z2) 
D3(0.z)  =  1 


(V-llb) 


Djd.z)  =  2-z 

03(2, z)  =  1(12  -6z  +  z2) 

D,(3,z)  =  1(120  -  60z  +  12 z^  -z^)  . 

j  6 


(V-llc) 
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VI.  SUMMARY 


The  effects  of  obscuration  on  target  acquisition  and  search  have  been  modeled 
mathematically.  In  Section  m,  a  statistical  mechanical  analog  between  the  alternating  states 
of  clear  and  blocked  lines  of  sight  on  the  one  hand  and  the  up  and  down  states  of  a 
ferromagnet  was  exploited  to  give  a  general  relationship  between  the  search  timelines  in  the 
presence  of  obscuration  and  the  free  energy  of  a  magnet  (Eq.  III-8).  Applying  this  general 
result  to  the  simplest  magnet  model  (the  Ising  model)  provides  an  explicit  solution 
(Eq.  rV-12).  This  choice  leads  to  exponential  distributions  of  the  obscured  and  clear 
intervals  (with  characteristic  exponents  of  and  p.*.,  respectively)  and  corresponds 
precisely  to  the  Monte  Carlo  model  of  obscurants  used  by  Kowalczyk  and  Rotman;  the 
model  parameters  are  determined  by  the  average  lengths  of  the  obscured  and  clear  intervals 
[(P-)~^  and  0+)~^]  and  the  unobscured  search  time  exponent,  a.  The  probability  of 
finding  a  target  increases  as; 

P«Kx*«i(0  =  1  -  e"^“'^“‘'^^‘'^[cos2  60  +  sin^  8ee‘^“‘]  (VI-l) 


where  p  is  the  overall  smoke  correlaticHi  inverse  time  constant  and 


cos^60 


=  (p  -  a<s>  )^  +  (1  -  <s>^) 


P  -  a  <s> 

fig 

2 


sin^66 

<s> 


p  -  a<s> 
Pa 

2 

P--  P. 


(VI-2) 


(<s>  is  the  "average  spin"  of  the  magnetic  analog;  the  average  visibility 
<‘n>  =  (1  +  <s>)/2.)  This  provides  explicit  expressions  for  the  numerical  results  of 
Kowalczyk  and  Rotman.  It  was  also  shown  that  this  model  could  be  equally  well 
represoited  as  a  3-state  Maikov  process. 
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In  Eq.  (VI- 1)  the  target  acquisition  task  was  modeled  as  a  cumulative  process  for 
which  the  total  time  available  on  target  was  the  relevant  quantity,  independent  of  how  that 
time  was  accumulated  or  how  many  times  the  line  of  sight  was  interrupted  by  obscurants. 
This  is  consonant  with  the  assumption  that  target  acquisition  is  the  consequence  of  a  series 
of  independent  glimpses,  each  of  which  is  of  shon  duration  compared  to  any  change  in  the 
visibility  conditions.  Thus,  the  task  is  performed  equally  well  when  10  seconds  of  clear 
viewing  is  provided  whether  it  is  provided  in  a  single  interval  of  10  seconds  or  in 
10  intervals  of  1  second  each.  Another  point  of  view  would  be  to  assume  that  a  single 
continuous  block  of  time  is  necessary  to  perform  a  task.  This  would  be  the  case,  for 
example,  for  a  command-guided  missile,  but  may  also  be  relevant  for  other  aspects  of  the 
target  recognidon  process.  The  Ising-Markov  model  that  was  used  to  model  the  cumulative 
process  can  also  be  used  to  model  the  probabilities  of  having  extended  (continuous) 
intervals  of  clear  line  of  sight  to  the  target. 

The  result  is  algebraically  complex  in  detail  [Eq.  (V-6)]  and  is  represented  by  an 
infinite  series  of  terms.  However,  numerically  the  results  are  relatively  simple  as  illustrated 
in  the  figures  (cf.  Figs.  V-l-V-3).  Fortunately,  a  simple  single  term  approximation  to  the 
exact  expression  can  be  derived.  The  probability  of  obtaining  an  interval  of  length  x  is 
given  by 


P+  =  1 


K(K-p) 

(K-P)(l+ . -) 

_ P-fP-  g-Kt 

2k  -  p  -  X  K  (K  -  P) 

(K-P)(l-^) 

_ P-t^  g-Kt 

2k-P-xk(k-P) 


(VI-3a) 


(VI-3b) 


where  P+  gives  the  probability  when  the  t  =  0  state  is  known  to  be  clear  and  P-  gives  the 
probability  when  the  initial  state  is  known  to  be  obscured. 
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The  parameter  k  is  determined  by  the  root  of  the  transcendental  equation^^: 

K(K-P)  +  p^P-e'^*‘*''^  =  0  .  (VI-4) 

This  single  term  approximation  is  valid  i<x  "large"  t/t  values;  in  practice,  the  error  2XXH-2 
is  negligible. 


There  are  two  real  roots  for  k;  the  remaining  (infinitely  many)  roots  are  complex.  One  root  of  this 
equation  is  always  Pf,;  the  root  to  be  employed  in  Eq.  (V-2)  is  the  other  real  root  In  the  q)ecial  case 
of  P+  -  P-  +  P+P-t  s  0,  the  real  roots  of  (VI-3)  are  degenerate  and  a  slightly  different  form  is 
use± 


P+ 


1  +  P^P-T^/2 


P- 


1  +  P^P_x2/2 
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